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In the vicinity of a phase transition, the order parameter starts fluctuating before vanishing at the 
critical point. The fluctuation regime, i.e. the way the ordered phase disappears, is a characteristics of 
a transition, and determines the universality class it belongs to. This is valid for thermal transitions, 
but also for zero temperature Quantum Phase Transitions (QPT) where a control parameter in the 
Hamiltonian drives the transition. 

In the case of superconductivity, the order parameter has an amplitude and a phase, which can 
both fluctuate according to well identifled scenarios. The Ginzburg-Landau theory and its extensions 
describe the fluctuating regime of regular metallic superconductors, and the associated dynamics of 
the pair amplitude and the phase. When the system is two-dimensional and/or very disordered, phase 
fluctuations dominate. 

In this article, we address the possibility that a new type of fluctuations occurs in superconductors 
with an anomalous dynamics. In particular we show that the superconducting to metal QPT that 
occurs upon changing the gate voltage in two-dimensional electron gases at LaAlOs /SrTiOs and 
LaTiOs /SrTiOs interfaces displays anomalous scaling properties, which can be explained by density 
driven superconducting critical fluctuations. 

More precisely, a Finite Size Scaling (FSS) analysis reveals that the product zu {u is the correlation 
length exponent and ^ the dynamical critical one) is zv ^ Zj2. We argue that critical superconducting 
fluctuations acquire an anomalous dynamics with z = 3, since they couple to density ones in the vicinity 
of a spontaneous electronic phase separation, and that = 1/2 corresponds to the mean-fleld value. 
This approach strongly departs from the conventional z = 1 scenario in disordered 2D systems based 
on long-range Coulomb interactions with dominant phase fluctuations. 

A zv = 2>l2 exponent has been recently reported in high Tc superconductor La 2 -xSrxGu 04 ultra-thin 
Aims upon gating, which could be explained by an analogous scenario. 

More generally, since superconducting fluctuations may couple to other electronic degrees of freedom 
such as nematic, charge-order, or density fluctuations, we anticipate that a concomitant abundance of 
soft critical modes will affect and characterize a whole class of two-dimensional superconductors. 
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The microscopic BCS theory of superconductivity and its phenomenological version, the Ginzburg-Landau theory 
in mean field approximation describe the emergence of a quantum coherent phase out of a regular Fermi liquid 
metal^, with a large electronic density, a well established Fermi surface and a nearly constant density of states at 
the Fermi level. Here the superconducting critical temperature Tc relates to the density of states at the Fermi level. 
Other superconductors such as high Tc cuprates or pnictides strongly depart from this picture, with non-Fermi liquid 
behaviors. In dilute superconductors such as doped SrTiOs Tc depends on the electronic density n, and not only 
on the density of states. 

In the traditional BCS-Ginzburg-Landau scheme the destruction of the ordered phase as the temperature approaches 
the critical temperature Tc is ruled by classic thermal fluctuations of both the amplitude and the phase of the order 
parameter. The spatial coherence length ^ in the fluctuating phase scales with the distance to the critical temperature 
with a exponent v oc {T — Tc)~^), while the dynamics of the fluctuations slows down with an exponent z = 2 
(r oc r being the characteristic lifetime of the fluctuations. This leads for instance to the standard Aslamasov- 
Larkin^ corrections to the conductivity above Tc. The same behavior is observed for a zero temperature Quantum 
Phase Transition (QPT), when a parameter K of the Hamiltonian is continuously changed so that superconductivity is 
destroyed, and the metallic phase restored beyond a critical value i^c^- In dirty two dimensional superconductors, the 
reduced screening of the Goulomb electron-electron repulsion may lead to pair breaking^ or to the direct localization 
of the superconducting pairs with Goulomb interactions rendering the superfluid density weaker and the conjugated 
phase fluctuations correspondingly softer. Thus, phase fluctuations dominate the physics of the system, which can be 
described by a spin XY model. In this context, Fisher argued that the QPT from a superconductor to an insulator is 
ruled by a dynamical exponent z = since long range Goulomb interactions are restored in these strongly localizing 
systems^. More recently, it has also been suggested that approaching the transition, disorder itself destroys coherence 
and inhomogeneities develop even in structurally homogeneous thin films. An intrinsic granular behavior emerges, as 
evidenced by local tunneling spectroscopy^’^, and explained by theories^’One can wonder how universal are these 
above presented descriptions of phase transitions in superconductors, with dynamical exponents z = 1 or z = 2. 

The recent discovery of a two-dimensional electron gas (2DEG) at oxide interfaces'^ opened fascinating perspectives 
since, in addition to strong confinement, they display d-electrons related electronic properties'^, a sizable Rashba type 
Spin Orbit Goupling^^ and quantum orders such as superconductivity^^ and magnetism^^. Quantum confinement of 
anisotropic d-orbitals leads to a rich and complex band structure^^’^^, which can be filled by electrostatic back-gating. 
A key point is that this band structure evolves with the filling through self-consistent shaping of the potential well^^. 
It was shown that, in some circumstances, the system displays a negative inverse electronic compressibility d[ijdn < 0 
(/i is the chemical potential), which leads to a spontaneous phase separation between high and low electronic density 
regions^^. The Rashba Spin Orbit Goupling can also foster an electronic instability and a correlated phase separation, 
which persists as the gate voltage is changed^^. Electronic correlations related to d-orbitals may also lead to negative 
compressibility as recently reported in LaAlOs /SrTiOs interfaces'^ together with phase separation^^’^^. A phase 
diagram emerges with an electronic phase separation region ended by a critical point. In this article, we show that 
a new kind of superconducting fluctuations takes place in the vicinity of this density critical point, with a z = 3 
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dynamics. 

Superconductivity at LaAlOs /SrTiOs and LaTiOs /SrTiOs epitaxial interfaces has been recently investigated. 
The analysis of the magnetic critical fields evidences the strong two-dimensional character of the superconducting 
2DEG which has a typical thickness th smaller than 10 nm, in agreement with self-consistent calculations of the 
electronic states in the quantum well formed at the interface^^, and a superconducting coherence length of ^ ^ 
40 — 70 nm {th < <f). With a Fermi energy in the 100 — 150 meV range^^’^^, the Fermi wavelength Ai? is ^ 10 nm^ that 
is slightly greater than the extension of the well (A^ > th). This 2DEG is therefore an extreme 2D superconductor, 
with an electronic sub-band structure due to the confinement. As the doping is changed electrostatically with a 
gate deposited at the rear side of the substrate (back gating), the superconducting temperature Tc is modulated, 
and shows a dome-like dependence with the gate voltage Vq ending with a Quantum Gritical Point. As Vq 

is increased, electronic sub-bands with increasing energy and spatial extension are populated. We showed that 
superconductivity is related to the occurrence of Highly Mobile Garriers setting at the edge of the quantum well, 
which correspond to high energy electronic states extending farther in the SrTiOs substrate^^. In agreement with 
direct superfluid density measurements on LaAlOs /SrTiOs interfaces'^, it is found that a small fraction (of the order 
of 10%) of the total electron density forms the superconducting state, thereby also supporting the scenario of an 
inhomogeneous superconducting state^^’^^. We recently studied the superconductor-localizing metal QPT that occurs 
when a perpendicular magnetic field is applied on the 2DEG^^. Through a finite size scaling analysis, we evidenced 
that the system is formed of superconducting puddles coupled through a metallic phase. Depending on the strength 
of the coupling controlled by Vq and the temperature, the critical behavior is dominated by the intra- or inter-puddles 
phase fluctuations of the superconducting order parameter. A striking result is that the exponent zn ^ 3/2 at the 
lowest temperatures, which departs from usual results in the literature (see Goldman’s review^^). However, recent 
measurements by Bollinger et al.^^ on a ultra thin high Tc superconductor film (1 to 2 unit cells thick), also found 
zu = 3/2 upon gating. 

In this article, we analyze the QPT of the superconducting 2DEG at LaAlOs /SrTiOs and LaTiOs /SrTiOs 
interfaces when using Vq as the control parameter of the Hamiltonian. The same zn ^ 3/2 exponent is found, which 
we interpret as a new kind of superconducting fluctuations controlled by an electronic phase separation critical point. 

Samples are grown by pulsed laser deposition of 15 unit cells of LaTiOs or LaAlOs on Ti02 -terminated SrTiOs 
substrate on the (001) direction (see Biscaras et al.^^, Herranz et al.^^ and Rastogi et al.^^ for details). A metallic 
back gate is evaporated at the rear of the SrTiOs substrate (500 jam thick) and connected to a voltage source (Vq )• 
Standard four probes resistance measurements are made with low current (10 nA) and low frequency (13 Hz) lock-in 
voltage detection. 

In Figure [1^, the sheet resistance Rs of a LaAlOs /SrTiOs sample is plotted as a function of temperature T 
for different gate voltages Vg -The superconducting state is progressively destroyed in favor of a weakly localizing 
metal^^’^^, as seen in the phase diagram Figure [lb. At the lowest temperatures, a plateau develops for a critical value 
of the sheet resistance Rc ^ 2.57 kQ/\3, which separates the superconducting and the normal phases. Plotting Rg as a 
function of Vq for different temperatures (Figure [2b) reveals a clear crossing point {Rc — 2.57 /cD/D, Vqc = —78.5 V) 
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which is a sign of possible QPT. A finite size scaling analysis (FSS)^ of the data in this region is made on Figure [2]3. 
All the sheet resistance data from 35 to 110 mK collapse onto a single function 


— = F 
Rc 


|Vg - Ksc 

2 ^ 1 /zv 


( 1 ) 


if zv ^ 1.6 ± 0.1. This critical exponent can then be retrieved (Figure [2t) by a scaling procedure^^ described in the 
methods section. 



Figure 1: LaAlOa /SrTiOs interface : (a) Resistance per square Rs as a function of temperature for different gate voltages Vb 
from +200F to — 200F. (b) Zoom in the low temperature region, where a plateau develops for Rc ~ 2.57A;fl/n. (c) Resistance 
per square Rs (left axis) and superconducting critical temperature Tc (right axis) as a function of the gate voltage Vb • The 
superconducting region (SC) is colored in blue. 


The behavior of LaTiOs /SrTiOs 2DEG is very similar. Figure [3^ displays Rg as a function of temperature T for 
different gate voltages : a critical sheet resistance Rc ^ 2.35 /Al/D separates the superconducting phase from the 
metallic one, and a crossing point is evidenced (Figure [Sb). The FSS analysis of the low temperature data points 
toward a QPT with a critical exponent zv ^ 1.6 ± 0.1 (Figure [3t). 

The exponent zv ^ 3/2 is somewhat surprising, and rarely found in the literature. The most common reported 
values are close to = 2/3, zv = A/oi zv = 7/3^^. Assuming a dynamical exponent z = the spatial exponent 
z/ = 2/3 would therefore correspond to a clean (2+l)D XY model, coherent with a scenario where phase fluctuations 
dominate the behavior of the system in the vicinity of the transition, with an enhanced dimension for the quantum 
character of the transition. Under the same assumption 2 : = 1, z/ = 4/3 and z/ = 7/3 are typical of a classical, 
respectively quantum percolating behavior in highly disordered systems. The main assumption in these arguments 
is that the dynamical exponent 2 ; equals 1, because of enhanced long range Coulomb interactions when the system 
becomes insulating^. Yazdani et al.^^ measured it in amorphous MoGe films and Markovic et al.^^ in a-Bi layers, and 
confirmed this value. 
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Figure 2: LaAlOa /SrTiOa interface : (a) Rs as a function of the gate voltage Vb for different temperatures from 0.035 to 
0.11 K. The crossing point is {Rc ~ 2.57A;fl/n,VGc = —78.51/). (b) Finite size scaling plot Rs/Rc as a function of |Vg — VgcI^ 
(see text for the dehnition of t). (c) Temperature behavior of the scaling parameter t (see text). The power law ht gives 
zu - 1.6 ± 0.1 

Following the same assumption ^ = 1, we were not able to perform a satisfying FSS with the above quoted u 
exponents in none of our samples. Moreover, since the non-super conducting phase is a weakly localizing metal and not 
an insulator^^’^^, screening effects should be substantial, and this might cast some doubts about the presence of long- 
range Coulomb forces customarily invoked to justify the z = 1 critical exponent. Indeed the highly metallic character 
of these materials leaves the possibility open that the z = 1 critical behavior usually invoked for the superconductor to 
insulator transition could be replaced by an over damped dynamics with z > 1. Thus, in the following, we investigate 
an alternative scenario, in which the dynamical critical behavior of the superconducting fluctuations is imposed by the 
coupling to nearly critical density fluctuations with z = 3. In this scheme, the experimental observation of zz/ ^ 3/2 
arises from a z = 3 dynamical critical behavior of the superconducting fluctuations together with a mean-field like 
exponent z/ = 1/2^^. 

Assuming that a phase separation takes place in the system as already mentioned^we sketch in Figure Ht 
the phase diagram of the electronic density as a function of the gate voltage. At low density (very negative Vq ), 
the electronic system is homogeneous, and becomes phase separated when entering the instability dome. In this 
region, static clusters of high density n 2 s are embedded in low-density regions of density ni^, whose proportions are 
given by the Maxwell construction. This intrinsic electronic inhomogeneity accounts for observed inhomogeneous 
superconducting properties^^’^^: indeed, if n 2 s is high enough, the islands can be superconducting^^, embedded in 
normal zone of density ni^. The dome ends at a Quantum Critical Point (QCP)^^, which is the analogous of the 
critical point of the classical liquid-gas phase diagram. In the vicinity of the QCP, critical density fluctuations are ruled 
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Figure 3: LaTiOs /SrTiOs interface : (a) Resistance per square Rs as a function of temperature between 0.03 and 0.2 K for 
different gate voltages from +301/ to —20V. A plateau is observed for Rc ~ 2.35kQ/n (b) Rs as a function of the gate voltage 
Vg for different temperatures from 0.035 to 0.15 K. (b) The crossing point is {Rc 2.35A;Q/n,VGc = —151/). (c) Finite size 
scaling plot Rg/Rc as a function of IVg — VgcI^ (see text for the definition of t). Inset : temperature behavior of the scaling 
parameter t (see text). The power law fit gives zz/ ~ 1.6 ± 0.1 

by a z = 3 dynamics. This substantially increases the effective dimensionality D^z and the system therefore displays 
a mean-field exponent u = 1/2^^. Let us first make the hypothesis that superconductivity is destroyed at the critical 
gate voltage Vgc of fho QCP . In other words, this voltage corresponds to the beginning of the phase separation and 
the occurrence of superconductivity (in high density droplets). In that case, superconducting fluctuations develop on 
density fluctuating droplets, and therefore acquire their dynamics with z = 3^^: consequently = 1/2, and zu = 3/2. 

There are good reasons for the density QCP to coincide with the occurrence of superconductivity. The first one is 
closely related to the physics of this system. The phase separation that takes place is likely due to non-rigid d-wave 
bands and sub-bands in the quantum well. The complex band structure at the interface is still under debate, but 
a consensus has emerged that the shape of the well, and therefore the spacing between sub-bands together with the 
energy of the bottom of the conduction band change with the filling of the well, and therefore the electron density. 
The more electrons are added, the deeper in energy is the well. This is the origin of the non rigidity of the bands. 
The lowest bands are of d^y character, and lie close to the interface. Upon doping, the populated upper band can be 
of dxy character^^ or dxztdyz one^^ or even a mixture between them through Spin Orbit Coupling^^. Whatever the 
scenario is, the common feature is that this band strongly delocalizes within the SrTiOs substrate when populated, 
which favors superconductivity^^’^^. In the framework of Ref. 21, the extension of the 2DEG as a function of gate 
voltage has been calculated (see Figure HJi). When entering the phase separated domain, the spatial extension of the 
band related to Usi does not evolve significantly while the one related to ns 2 increases rapidly beyond 4 nm, where 
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Figure 4: (c) : Sketch of the electronic phase separation at oxide interfaces. A phase separation line ending by a density QCP 
marks out the homogeneous from the phase separated zone (in purple) in the Vb vs Ug phase diagram. For a sample with a 
density Ug (red arrow) at Vg < Vgc, the system has an homogeneous density. When further increasing Vg , it enters the phase 
separated zone, and decomposes into low density Ugi (yellow broken line) and high density ng 2 (blue broken line) droplets, as 
sketched within the dome. In the low density droplets, carriers are localized near the interface and remain normal (a), while 
in the high density ones, carriers extend far from the interface (b) and display superconductivity. The potential well presented 
in (a) and (b) have been calculated using Biscaras et al. method described in ref 20, for gate voltages of - 200 and + 200 V 
respectively, (d) : Within the Scopigno et al. theory^^, calculated gas extension d as a function of Vg • Beyond Vgc (which 
for the parameters chosen in the theoretical calculations occurs above Vgc ~ — lOOV), d remains constant for the low density 
regions at n = Ugi (yellow dots) , and increases in high density droplets at n = ng 2 (blue dots). 

the superconducting phase shows up^^. Consequently, phase separation and superconductivity appear in the same 
region of the phase diagram. In addition, superconductivity can, by itself, favor phase separation when Tc depends on 
the electronic density as it is the case here. If the pairing energy A is a function of the system may gain energy 
by forming superconducting regions of high density in non superconducting low density ones. A Ginzburg-Landau 
approach extended to a diffusion model confirms this possibility, and the occurrence of an electronic instability in this 

OQ 

case . 

Conversely, the electronic compressibility is strongly enhanced near a density-driven superconducting critical line 
Tc(n)^^. Thus, the QCP related to the phase separation dome (where the compressibility diverges), is attracted by 
the superconducting QCP (where Tc{n) vanishes), and the two QCP tend to merge. 

Seldom have been reports of 2D superconductor to weakly localized metal or insulator transitions controlled by the 
electronic density so far. Focusing on critical exponents extracted from a complete FSS analysis (see Supplementary 
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material), a few results are found in the literature. Parendo et al.^^ reported zv ^ 2/3 in amorphous Bi films, 
corresponding to a clean (2+l)D XY model in this highly metallic conventional s-p compound. In 7 u.c. thick 
YBa 2 Cu 307 layers^^, zv ^ 2.2 has been found corresponding to quantum percolation, which is not surprising in highly 
disordered thin films with up to 6 u.c. of dead layer. However, in 1 to 2 u.c. thick La 2 -xSrxCu 04 samples, Bollinger et 
al. found exactly zv = 3/2^^, while Garcia-Barriocanal et al. reported zv from 1.4 to 1.8 for 4 u.c La 2 -xSrxCu 04 thin 
films^^, not far from it. This system shares similarities with the LaAlOs or LaTiOs /SrTiOs interfaces, such as an 
extreme 2D character (thickness of a few unit cells) and a low electronic density at the transition (in the 10^^ e~ jcnr? 
range). It is worthwhile mentioning that the magnetic field driven transition looks also very similar in the two systems, 
with two critical regimes in temperature^^’^^. There is no report of negative compressibility induced phase separation 
in La 2 -xSrxCu 04 , but other intrinsic electronic inhomogeneities have been found such as fluctuating Charge Density 
Waves The coupling of the superconducting fluctuations to other critical electronic modes near Charge Density 
Waves^^, stripes or fermionic nematic phases^^ QCP may lead to anomalous dynamics as well. It would be interesting 
to explore more widely this possibility in cuprates, Fe-pnictides or other exotic superconductors. 

In conclusion, we studied the superconductor to metallic quantum phase transition in LaAlOs /SrTiOs and LaTiOs 
/SrTiOs interface as a function of the electronic density tuned by a gate voltage. The critical exponents product 
zz/ ^ 3/2 is compatible with density driven fluctuations, where the superconducting fluctuations are coupled to density 
ones in the vicinity of an electronic phase separation critical point. This new type of fluctuations may be observed 
in other 2D superconductors, and perhaps also in superfluids such as on aerogels, where a similar zv has been 
found^^. 
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Methods 

The resistance is rewritten as R{6,t) = RcF{\6\t) with t an unknown parameter that depends only on T, and 
^ = (Eg — Vgc), the distance to the critical point. The parameter t is then found at each temperature T by optimizing 
the collapse around the critical point between the curve R{SR{T)) at temperature T and the curve R{SR{To)) at the 
lowest temperature considered Tq, with t(To)=l. The dependence of t with temperature should be a power law of the 
form t = in order to have a physical sense, thus giving the critical exponent product zv . The interest of 

this procedure is to perform the scaling without knowing the critical exponents beforehand. 
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